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Abstract 

The construction of a nonautonomous mixed mKdV/sine-Gordon model is proposed by 
employing an infinite dimensional affine Lie algebraic structure within the zero curvature 
representation. A systematic construction of soliton solutions is provided by an adaptation 
^ of the dressing method which takes into account arbitrary time dependent functions. A 

particular choice of those arbitrary functions provides an interesting solution describing 
\^ the transition of a pure mKdV system into a pure sine-Gordon soliton. 
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1 Introduction 



Sometime ago, the study of nonlinear effects in lattice dynamics under the influence of a weak 



CSJ dislocation potential has lead to a mixed mKdV/sine-Gordon equation pj. The system was 
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QQ shown to admit multisoliton solutions and an infinite set of conservation laws [Ij. More recently 

^ the two-breather solution was discussed in p] in connection with the propagation of few cycle 
O 

O pulses (FCP) in non linear optical media. According to ref. [2] the general mKdV/sine-Gordon 



equation, in fact, describes the propagation of a ultrashort optical pulses in a Kerr media. 
Moreover, it was shown in [3] that, when the ressonance frequency of atoms in the physical 
system are well above or well below the characteristic duration of the pulse, the propagation is 
described by the mKdV or sine-Gordon equations respectively. The main object of this paper 
is to provide a systematic construction of soliton solutions that describe the transition between 
the two regimes, i.e. governed by the mKdV and sine-Gordon equations. This is accomplished 
by considering the mixed integrable model proposed in [1] with two arbitrary time-dependent 



coefficients. In this paper we show the integrabihty of the mixed model with time dependent 
coefficients and that, by suitable choice of these coefficients as a smooth step-type functions 
(as shown in figs. 1 and 3) we obtain exact solutions for the mKdV-SG transition and hence a 
more realistic description of such phenomena. 

2 Algebraic Formalism 

In ref. the algebraic structure of the mixed mKdV/sine- Gordon equation was formulated 
within the zero curvature representation and a graded infinite dimensional Lie algebraic struc- 
ture as we shall now briefly review. Consider the associated Q = sl{2) Lie algebra with gen- 
erators satisfying [h,E±a] = ^'^E±a, [Ea,E_a] = h and grading operator Q = 2A^ + ^h. Q 
decomposes the associated affine Lie algebra sl{2) into graded subspaces, Q = (BiGi, 

G2m = {X'^h}, = {4'^"^'^ = X^ (E^ + XE_^) , E^-"^^'^ ^ A™ {E^ - J}, (1) 

m = 0, ±1, ±2, .... In [1], a simple proof that a mixed mKdV/sine-Gordon hierarchy is indeed 
an integrable model follows from the zero curvature representation of the integrable hierarchy 
generated by 

[d, + E« + Ao, dt, + + + ■ ■ ■ + D^-'^ = 0, (2) 

where D^^^ E Qj and Aq = vh contains the field variable v = v {x,t). According to the subspace 
decomposition ([2]) for = 3 and t = which corresponds to the mixed mKdV-SG equation. 
Let us parametrize, 

D^'^ = a,E^^^ + hE^^\ D^'^ = C2XK 

= a_i4-') + 6_iEi-'). (3) 
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The grade by grade decomposing of eqn. g leads to 

63 = 0, d^as = 0, C2 = a^iV, bi = ^d^C2, (4) 
d^ai + 2vbi = 0, + 2vai - 2co = 0, (5) 

d^a_i + 2vb_i = 0, d^b_i + 2va_i = 0, (6) 

together with the equation of motion 

d.co - dtv - 26_i = 0. (7) 

In solving eqns. Q we find 

03 = aaW, h = ^^'^^v^, C2 = a3{t)v, (8) 
where a3(t) is an arbitrary function of t. Introducing ^ in the first eqn. ([s]), we obtain 

/ 



which implies that 



dx lai + a-iit)—] = 0, 



where /i(t) is another arbitrary function of t. It therefore follows that 

ai = /i(t)-a3(t)^. (9) 

Substituting ^ in the second eqn. ([s]), we get 

co = ^{vxx-2v')+h{t)v. (10) 

Adding and subtracting eqns. ([6]), we obtain 

dxa± = =F2t;a±, (11) 
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where we have denoted 

a± = a_i ± 6-1. 



Without loss of generahty we may solve (11 ) by changing the variable 



vh = -d^BB-^ = ^^h, B = e"^^, (12) 

which leads us to a± = f^i{t)e^'^'^, where f-i{t) is another arbitrary function of t. Writing 

a+ + a_ a+ — a_ 
«-i = — ?; ) = ?; ) 



we find 

a-i = f^iit) cosh (2</)) , 6_i = sinh (2</)) . (13) 

Substituting (10), (12) and (13) in ([T]), we finally obtain 
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- Q(g(t)xx) + /i(t)0xx - 0xt + 2/_i(t) sinh (20) = 0. (14) 



Considering /i(t) = 0, a-i{t) = constant and /_i(t) = constant we find the usual mixed 
mKdV/sine-Gordon equation. For fi{t) = 0, 03 (t) a given numerical constant 7^ we recover 
eqn. (10) of ref. [5j. Moreover for f-i(t) = 0, we recover equation considered in ^ with a 
choice of coefficients that makes the model integrable. 

We should point out that by change of coordinates (see for instance |7j) (x,t) — )■ {x,t) = 
{x+V(t),t) where Vt = fi{t) followed by a subsequently change i ^ T = J a^{t)dt and re-scaling 
/_i — 7- r/ leads to 

\{4>xxxx - Qfyxx) - <Pxt + 2f/(t) sinh (20) = 0. (15) 



Although eqn. (15) corresponds to the equation discussed in [5] the object of this paper is to 



consider a class of solutions that interpolates between the mKdV and sine-Gordon equations. 



This is more conveniently accomplished by employing eqn. ( 14 ) where the two arbitrary func- 
tions asit) and f-i{t) (with fi{t) = 0) can be chosen as step-like limiting functions (Figs. 1 
and 3) as we shall see. 

3 Construction of Soliton Solutions 

In order to construct, in a systematic manner, the soliton solutions of the mixed model let us 
now recall some basic aspects of the dressing method (see for instance p]). The zero curvature 
representation ^ implies in a pure gauge configuration, i.e., 

d^ + E + Ao = d^TT-\ dt + /^^^^ + ■ ■ ■ D^-^^ = dtTT-\ (16) 

In particular, the vacuum is obtained by setting (j)yac = Q which implies, 

a.ToTo-i = -E'i\ dtToT,' = -a3(t)4'^ - /-i(t)4"'^ - /i(t)4'^- (17) 

which after integration yields 

To = exp (^- y* dt'a^{t')Ef - j' dt' f-i{t')E^+^^ - j' dt' h{t')E^+^^ exp(-x4^), (18) 

Following the dressing method explained in [8] and employied in we define the tau-functions 

Tn = {\n\B\\n) = {\n\TogT^'\Xn) , (19) 

where Xn,n = 0, 1 are fundamental weights of the full affine Kac-Moody algebra s/(2), g is 
a constant group element which classifies the soliton solutions and S is a zero grade group 



^For a general member of the hierarchy evolving according t = <2n+i, the vacuun configuration implies 

n 

fc=i 
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element containing the physical fields. In order to ensure heighest weight representations we 
now introduce central extensions within the affine Lie algebra, characterized by c, i.e., 



and define highest weight representations, i.e. 



/i|An) = 5n,l|A„), c|A„) = |A„), ^i|A„) = 0, 2 > 0, 



(20) 



n = 0, 1. Under this affine picture the group element B acquires a central term contribution. 



B = e-^^e 



cj)h —vc 



(21) 



In order to obtain explicit space-time dependence from the r.h.s. of (19) we consider the vertex 



operators. 



\/ (7) = E (^"^ - ^c5„,o)7-^" + ^(^"+1)7-2-1, 



satisfying 



4'"+'\v^(7) =-27^"+V(7). 



2n+li 



For a general M-soliton solution the group element g m (19) is written as 



M 

i=i 



where ctj are arbitrary constants. We therefore obtain 

M 

To = e-^ = (Aoine"^''^'^"'*^''^^^'^|Ao), 

M 

n = e"^-'^ = (All ne°^''''^"'*^^^^''^|Ai), 
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(22) 



(23) 



(24) 



where 

t) = e27.-+27|A3(i)+27,Fi(<)+2,-iA_,(t)^ (25) 



As{t) = I dta^it), Flit) = I dtfiit), A_i(t) = I dtf-iit). 



As an illustrative example, we consider the one and two soliton cases, M = 1, 2, where 



l—sol 

'o 



1— soZ 
1 



1 , "1 

1 + 



(26) 



and 



2-sol 







.2-soZ 



Oil OLi 

1 - yPi - + aia2^i,2PiP2, 
+ yPi + yP2 + ai«2^i,2PiP2, 



(27) 



respectively. In order to obtain (26) and (27) where we have used the fact that 



(A„|V^(7)|A„) = bn, 



(A„|V(7i)l^(72)|A„) = A 



1,2 



1 / 7i - 72 
4 Wi +72, 



n = 0, 1 

\ 2 



The general two soliton solution can then be written as 



In 



In 



1 - f-Pl - ^P2 + ai«2^1,2PlP2' 

1 + ^pi + + aia2^i,2piP2, 



(28) 



while the one soliton is obtained from (28) by setting (X2 = 0. 



^ In considering a general (2rt + l)-th member of the hierarchy, 
Pj{x,t) = e " ^ ;t z^^^i iwT ,j A2„+i(i) = / (ita2„+i(0, ^2fc-i 



dtf2k-i{t). 



7 



4 Applications 



According to ref. [2] the propagation of a FCP with frequency a; on a dielectric media with 
characteristic frequency Q, u << Q can be described by the mKdV equation 

where the coordinates z and r correspond respectively to the propagation distance and retarded 
time, while the electric field E = (f)^. For the case where u » fl, the system is described by 
the sine-Gordon equation, 

— 6sin0 = 0. 

If we now consider a dielectric media with two characteristic frequencies Qi, Q2, the case in the 
regime fii << a; << ^^2, is described by the mixed mKdV-SG equation , 

4>ZT + a Q^r^rr + <Ptttt^ - 6 Slu = 0, 

where the two constants a and b are related to the non-linear and dispersion properties of the 
media. 



In order to adapt model (14) to such situation, define 



re-scaling — > |0, t — > 2;, x — r, eqn. (14) becomes 

3 



aOliz) (^(l)rrrr + 1^4>UrT^ + " &^2(^) Slu = 0. (29) 



If we now substitute — )■ —2iak, and make use of the identity 



2i \l-iX 
8 



we find that the two sohton solution (28) may be written as 



where 



/ aipi + a2p2 \ 

4 arctan 

\l- 4aia2Ai,2PiP2 / 



Pj = exp 



(30) 



Asiz) = -4a J dz'Siiz'), 
4.1 Transition mKdV-SG 



A_Az) 



fdz'e^iz') 



b r 

4 



(31) 



Consider now a dielectric media in which fii >> 27^ in the region z < Zi and 1^2 << 27^ in 
the region z > Z2 with zi > Z2 such that there exist an overlap region in which the media 
admits the two characteristic frequencies, ^2 << '^Ij << ^^i- As an example to describe such 
a realistic situation, we take both Oi{z) and ^2 (-2) as step- like functions (see fig. 1) with 

6i{z) = - - arctan [Pi{z - Zi)] , e2{z) = + - arctan [P2{z - Z2)] , 

2 TC 2 TT 

where /3i e P2 are phenomenological parameters describing the transition between the two 
medias (see fig. 1). 



1 n 








0.6 




0.4 




0.2 






Figure 1: Plots of Oi^z) and ^2(2) for /3i = /32 = {f , f , 27r}, with z^ = 2 and Z2 = -2 
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They can therefore be integrated from ( 32 ) to yield 



4a Z 71 Znpi L 



-A-,{z) 



Z [Z - Z2 

2 



arctan [/32{z - Z2)] - 5-ln 
TT 27rp2 



1 + PI(Z-Z2) 



The f3i, (32 — )■ +00 hmit, for zi = Z2 = correspond to the system governed by the pure mKdV 
in the region z < and by the pure sine-Gordon equation for z > 0. In such hmit, we have, 



Mz) = \ 1 



\z\ 
z 



1 



(32) 



and 



^A_^{z) = ]^{z+\z\). 



(33) 



Figure 2 below shows the transition mKdV-SG for the one soliton solution, ai = 1, 0:2 = 
0. The plot on the right shows the soliton solution viewed from the above and displays the 
transition (different velocities) from the mKdV to the SG solitons. 





Figure 2: Plot (left) and contour plot (right) of (pr with a = ^, 6 = — |, [3i = [32 = Svr, zi = 1, 
Z2 = -1 and 7i = ^. 
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4.2 Transition mKdV-SG-mKdV 

Another example consist in two equal media separated by a second one describing, for instance 
the mKdV-SG-mKdV transition. Mathematically this situation may be described by combining 
theta-type functions (see fig. 3), i.e., 



9Az) 



1 arctan 

TT 



f3i{z - Zi) 



arctan 



TT 



(32{Z - Z2) 



Zl < Z2, 



02{z) = — arctan [Pi{z — zi)] arctan [(32{z — Z2)] , Zi < Z2, 

TC TT 

with Zl < Zl, and Z2 < Z2 to guarantee the existence of an overlap region between the two 
medias. 

B(z} 



Figure 3: Plot of Oi{z) (dashed) and 6*2 (-2) (continuum) with (3i = (32 = (3i = (32 
Z2 = 6, Zl = —4 and Z2 = 4. 



in, Zl 



After integration (32) we find 
1 



4a 



-A,(z) 



z — arctan 



TT 



(3i{z - Zl) 



H arctan 



TC 



(32{Z - Z2) 



1 



-A.i{z) 



2n(3i 

z - Zi) 



In 



l + Pl{z--Zif 



1 



In 



TT 
1 



arctan [(3i{z — Zi)] 



2ix(32 

- Z2] 
TT 



271(3, 



In 



1 + Pliz - Zif 



+ 



271(3, 



In 



i + i3i{z~z2y\, 

arctan [(32{z — Z2)] 
l + (3l{z-Z2) 



The limit (3i, (32, (3i, (32 — )■ +00, when zi = zi e Z2 = Z2 with zi < Z2, corresponds to the pure 
mKdV case in the region z < zi and z > Z2, and pure sine-Gordon in the region zi < z < Z2. 
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Under such limiting case we have 
9Az) = 1 



1 
2 



Z — Zi 




Z - Z2 




[z-zi 


) {^-^2) J 



9,{z) 



1 ( \z — Zi\ 



\Z - Z2\ 



2\{z- zi) {z - 2:2), 



such that 



^M^) = Z -^i\z - Zi\-\Z - Z2\) 



4 1 

-A.l{z) = -{\Z - Zi\-\Z - Z2\) 



Figure 4 represents the niKdV-SG-mKdV transition for one sohton solution. 





Figure 4: Plot (left) and contour plot (right) of with a = j^,b 
-4, Z2 = A, zi = -3, Z2 = 3 and 71 = j^. 



-|,/3i = /32=A = /92 = 87r, 



4.3 Two soliton solution 

The two soliton solution (ai = 02 = 1) can be represented by 

As a conclusion, we have adapted the general dressing construction of soliton solutions to 
the mixed mKdV-SG hierarchy with arbitrary "time" dependent functions. The choice of such 
arbitrary functions as step-type functions allowed exact solutions describing smooth transitions 
from the mKdV to sine-Gordon regime and therefore more realistic models. 
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Figure 5: Plot (left) and contour plot (right) of (f)r with a = T^,6 = — |,/3i = /32 = /3i = /32 = Svr, 



zi = zi = -4, 2:2 = 2:2 = 4, 7i = I e 72 - 
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